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ABSTRACT 

We present the WFC3 Infrared Spectroscopic Parallel (WISP) Survey. WISP is obtaining slitless, 
near-infrared grism spectroscopy of ~ 90 independent, high-latitude fields by observing in the pure 
parallel mode with Wide Field Camera-3 on the Hubble Space Telescope for a total of ~ 250 orbits. 
Spectra are obtained with the G102 (A = 0.8 — 1.17 /im, R~ 210) and G141 grisms (A = 1.11 — 1.67 /im, 
R ^130), together with direct imaging in the J- and H-bands (FllOW and F140W, respectively). In 
the present paper, we present the first results from 19 WISP fields, covering approximately 63 arcmin^. 
For typical exposure times (~ 6400 s in G102 and ~ 2700 s in G141), we reach 5cr detection limits for 
emission lines of / ~ 5 x 10~^^ ergs s~^ cm~^ for compact objects. Typical direct imaging 5<j-limits 
are 26.8 and 25.0 magnitudes (AB) in FllOW and F140W, respectively. Restricting ourselves to the 
lines measured with highest confidence, we present a list of 328 emission lines, in 229 objects, in a 
redshift range 0.3 < 2: < 3. The single-line emitters are likely to be a mix of Ha and [GUI] 5007,4959 
A, with Hof predominating. The overall surface density of high-confidence emission-line objects in our 
sample is approximately 4 per arcmin^. These first fields show high equivalent width sources, AGN, 
and post starburst galaxies. The median observed star formation rate of our Ha selected sample is 
4 M0 year~^. At intermediate redshifts, we detect emission lines in galaxies as faint as J^_140 ^ 25, 
or Mr < — 19, and are sensitive to star formation rates down to less than 1 M© year~^. The slitless 
grisms on WFC3 provide a unique opportunity to study the spectral properties of galaxies much 
fainter than L* at the peak of the galaxy assembly epoch. 
Subject headings: galaxies: distances and redshifts - galaxies: statistics - infrared: galaxies - surveys 



1. INTRODUCTION 

While there is considerable uncertainty regarding the 
detailed star formation history of galaxies of various 
masses and types, it has become clear that the span 
of cosmic time 0.5 < z < 2.5 encompasses most of 
the star formation in the history of the Universe (e.g. 
iHopkins. A.1l2004l : iDaddi et al.l[ 2QQ7^). This broad epoch 
is also when many galaxies were "assembled," when ma- 
jor mergers and the accretion of satellite galaxies and gas 
combined to produce the variety of morphologies and 
structures of modern galaxies. There is mounting ev- 
idence for an evolution of star-forming galaxies during 
this time period, from high-mass galaxies dominating at 
2; ~ 3, to low- mass galaxies at the present. Deep, wide 
IR surveys are required to cover much of the peak epoch 
of star formation, because of the redshifting of the key 
diagnostic spectral features of stellar evolution, and to 
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minimize the effect of obscuration by dust. 

The difficulty of obtaining observed-frame near- 
infrared spectra causes much of our empirical perspec- 
tive on galaxy evolution to be biased by the way in 
which spectroscopic samples are selected. Color selec- 
tion reliably selects common objects, but misses unex- 
pected ones; magnitude-limited samples necessarily fa- 
vor bright-continuum objects. Ground-based searches 
for emission lines from faint high-redshift galaxies are 
severely impacted by the bright NIR airglow. Ex- 
cluding OH terrestrial airglow frequencies restricts red- 
shift coverage and compromises redshift identification. 
Space-based spectroscopic surveys offer significant ad- 
vantages for unbiased target selection. The grisms on 
the Near Infrared Camera and Multi Object Spectrom- 
eter (NICMOS), the Space Telescope Imager and Spec- 
trograph (STIS) and the Advanced Camera for Survey 
(ACS) have demonstrated the power of spectroscopy 
from above the atmosphere, wher e the low background 
emission permits slitless obs erving (iGardner et al. 1998[ 
McCa rthv et al. 1119 99: Tep litz et al.l 120031 : iPirzkal et al.l 
2004; Drozdovsky et al.. .20051 ). Furthermore, the ter- 
restrial wavelength gaps can only be filled from space, 
as first demonstrated by the pure-parallels program 
with the NICMOS G141 grism (McCarthv et al. 1999; 
Yan e t al. II1999I : iShim et al.l l2009). 

Much as ground-based objective prism surveys (e.g. 
Gallego et al. 1999) were used to define fiux-limited sam- 
ples of Ha emitters at z ~ 0, the NICMOS grism observa- 
tions produced complete samples of emission-Hne objects 
at intermediate redshifts. Yan et al. (1999) and Hopkins 
et al. (2000) used observations with the NICMOS H- 
band (G141) grism to derive Ha luminosity functions. 



and to show that the luminosity density, and hence star 
formation rate density, has evolved by a factor of ~ 10 
from z ~ 1.4 to the present. Shim et al. (2009) used 
a larger sample to derive improved Ha luminosity func- 
tions in two redshift bins and confirmed the trends seen 
in the analysis by Yan et al. and Hopkins et al. The 
complete reductions/extractions of grism data for the 
principal modes of NIC MOS were released in the Hubble 
Legacy Archive fHLA. iFreudling et al.ll2008[ ). Nearly ah 
of the published results from the NICMOS grisms were 
based on data with the G141 H-band disperser. The poor 
sentivity of the NICMOS detector at short wavelengths 
limited the utility of the J_110-band grism. 

The grisms on ACS have been used to discover faint 
Lya emission sources and faint Lyman continuum break 
sources at 2; ^ 3 — 4 ( Malhotra et al. 2005 ; Pirzkal et al. 
l2007f ) as well as to derive spectroph otometric redshifts for 
a large sample of faint galaxies (e.g. lStraughn et ani2008l 
[2009). Extremely deep ACS grism observations have 
shown the power of low resolution spectroscopy above the 
atmosphere to detect spect ral breaks indicative of pas- 
sively evolving svstem s fe.g. iDaddi et a l. 2005) or active 
star f orming galaxies (*Pirzkal et al.ll2Q04i: .Rhoads et al.. 
120051 : iMalhotra etld] 2005; P asauali et al.l I2006D . The 
full release of the reduced ACS grism data in the HLA 
is underway (Kuemmel et al. 2009b). The STIS slit- 
less modes work at higher resolution than the NIC- 
MOS and ACS grisms, but the relatively low sensitiv- 
ity of STIS has limited its ability to sample populations 
beyond the reach of ground-base spectroscopic surveys 
(jTeDlitz et al. II2003D . 

Wide Field Camera Three (WFC3) on Hubble has a 
larger field of view, significantly better detectors and su- 
perior sampling compared to camera 3 on NICMOS. The 
net improvement in survey speed is approximately a fac- 
tor of 20. This enhanced survey speed does not fully 
capture the gains provided by WFC3. The quality of 
the spectra, in terms of resolution, a well-focused point 
spread function, and uniform sensitivity are also greatly 
superior to those offered by NICMOS. 

In HST's Cycle 17, we began a new observing pro- 
gram of WFC3 pure parallels, the WFC3 Infrared Spec- 
troscopic Parallel (WISP) survey. In ~ 250 orbits, the 
program will obtain near-infrared slitless spectroscopy of 
~ 90 uncorrelated high latitude fields (PI=Malkan, GO- 
11696). This will produce a large sample of emission-line 
galaxies to measure the star formation rate continuously 
from at z ~ 0.5 to z ~ 2.5, over which ground-based 
searches are severely limited. The study of dust and 
metallicity in a large, unbiased sample of galaxies and 
the evolution of the star formation density are among 
the main goals of our project. In this paper, we present 
the analysis of the first 19 fields, from which we have se- 
lected high-confidence emission lines. In Sections [2] and 
[3l we describe the survey and data reduction. We discuss 
the results in Section HI and finally offer some thoughts 
on the future possibilities in Section [5l Throughout, we 
assume a A-dominated flat universe, with J^o = 71 km 
s~^ Mpc~^, l^A = 0.73, and ^^n = 0.27. All magnitudes 
are in AB system. 

2. OBSERVATIONS 

The infrared channel of the WFC3 utilizes a IKx IK 
HgCdTe detector array. At a plate scale of 0.13 arcsec 



px~^, the total field of view covered is 123'' x 136''. There 
are two IR grisms on WFC3: G102 covering the 0.8 — 1.1 
fim range with a dispersion of 0.0024 /im/pixel and G141 
covering the 1.07 — 1.7 //m range with a dispersion of 
0.0046 /im/pixel. These correspond to average resolving 
powers for point sources of i? = 210 for the G102 and R = 
130 for the G141. For an unresolved line, the FWHM 
will be about 2 pixels. The detai ls of the grisms and 
their characteristics can be found in lKimble et al.l (|2008f ) 
and on the WFC3 website^. 

Each target field was observed briefly in direct imaging 
mode, using the broadband filter that matches the grism 
spectral coverage most closely: FllOW for the G102, and 
F140W for the G141. As detailed in the next section, 
all of our observations were carried out in the "pure- 
parallel" mode with the Cosmic Origins Spectrograph 
(COS) or STIS operating as the prime instrument. 

2.1. Pure Parallel Observation Scheduling 

Pure-parallel observing with HST works as follows: 
Observing programs approved on COS and STIS con- 
stitute a parallel-observing 'opportunity' (see Tabled]). 
Observations with these two instruments, by necessity, 
require stable, long integrations, which means that the 
entire focal plane of HST is held at a single pointing and 
fixed orientation angle for long periods of time. Thus 
it is simply a matter of turning on one of the other in- 
struments in the focal plane to obtain quality data in 
a parallel field, offset by a 5.5' and 4.75' from the COS 
and STIS primary target, respectively. This motivated 
the strategy of a number of programs using WFC3 in par- 
allel. Two other HST Cycle 17 programs were approved, 
both to do infrared direct imaging (GO- 11702, PI = Yan; 
GO- 11 700 PI = Trenti). Pure-parallel observing is per- 
fectly suited to the goals of the WISP Survey, since it 
allows us to obtain deep, continuous, IR spectroscopy in 
dozens of uncorrelated fields. 

We selected our parallel targets with a preference for 
long integration times and high-galactic-latitude (> 30 
degrees out of the galactic plane) fields. In order to 
reach similar depths in both the G102 and G141 spec- 
tra, we planned an ~2.5:1 G102:G141 integration time 
ratio. Additionally, the process of extracting the slitless 
spectra from the grism images requires direct imaging 
of the same fields (see Section [3]) and so we obtained 
imaging in FllOW and F140W (corresponding to J_110- 
and H_140-band, respectively) on the same orbits as the 
grism imaging with an ~6:1 grism:direct imaging inte- 
gration time ratio. To date, we have received, reduced, 
and analyzed 19 fields (Table [I]), the results of which are 
presented here. 

Cycle 17 was the first year a powerful new procedure 
was adopted for scheduling and optimizing parallel ob- 
servations. Unlike previous cycles, each of the Paral- 
lels teams was allowed to select individual observing op- 
portunities, with accurately known durations of each or- 
bital visit, well in advance. Since every parallel visit is 
unique, with its own combination of full- and partial or- 
bits, we optimized each individual orbit to obtain the 
exposure time ratios mentioned above in the available 
time, thereby reaching a carefully balanced set of sensi- 
tivities. This tremendous improvement over earlier years 

^ http://www.stsci.edu/hst/wfc3/ 



TABLE 1 

Schedule of Observations to date for the WISP Survey 



Field 


Date 


RA 


DEC 


Number of 


FllOW 


G102 


F140W 


G141 


Primary 


Primary 


Primary 






[HMS] 


[DMS] 


Orbits 


[Sec] 


[Sec] 


[Sec] 


[Sec] 


Proposal ID 


Visit No. 


Instrument 


1 


2009/11/24 


01 06 35.29 


+15 08 53.8 


4 


884 


4815 


506 


2609 


11720 


01 


COS 


5 


2009/12/20 


14 27 06.64 


+57 51 36.2 


5 


1034 


5515 


1034 


5515 


11720 


13 


COS 


6 


2009/12/24 


01 50 17.18 


+13 04 12.8 


4 


609 


3609 


862 


5015 


11727 


09 


COS 


7 


2009/12/25 


14 27 05.91 


+57 53 33.7 


7 


834 


6318 


1112 


6224 


11720 


14 


COS 


9 


2009/12/29 


12 29 44.31 


+07 48 23.5 


4 


759 


4612 


684 


3712 


11703 


01 


STIS 


10 


2010/01/02 


09 25 07.84 


+48 57 03.0 


3 


631 


3909 


406 


2209 


11742 


06 


COS 


11 


2010/01/03 


11 02 17.38 


+10 54 25.4 


3 


556 


3709 


456 


2006 


11742 


07 


COS 


12 


2010/01/04 


12 09 25.25 


+45 43 19.8 


5 


1312 


8221 


606 


3009 


11741 


04 


COS 


13 


2010/01/09 


01 06 38.77 


+15 08 26.2 


3 


556 


3009 


506 


2409 


11720 


03 


COS 


14 


2010/01/09 


02 34 56.80 


-04 06 54.5 


4 


834 


6215 


481 


2809 


11741 


03 


COS 


15 


2010/02/12 


14 09 42.47 


+26 21 56.0 


5 


1612 


8321 


531 


2609 


11741 


26 


COS 


16 


2010/02/16 


02 34 54.72 


-04 06 42.5 


5 


1087 


6921 


584 


2509 


11741 


01 


COS 


17 


2010/02/18 


02 13 38.11 


+12 54 59.3 


3 


534 


3409 


559 


3409 


11727 


11 


COS 


18 


2010/02/19 


12 29 17.25 


+10 44 00.6 


5 


534 


3512 


762 


3824 


11561 


04 


COS 


19 


2010/02/20 


02 34 54.29 


-04 06 30.5 


5 


1187 


8721 


484 


2809 


11741 


02 


COS 


20 


2010/02/21 


14 09 41.15 


+26 22 15.1 


7 


1815 


8430 


559 


2812 


11741 


25 


COS 


21 


2010/03/05 


09 27 55.77 


+60 27 05.3 


1 








353 


2006 


11728 


10 


COS 


22 


2010/03/05 


08 52 44.99 


+03 09 09.6 


1 








253 


1806 


11728 


13 


COS 


23 


2010/03/06 


09 43 16.12 


+05 27 37.1 


3 








681 


4115 


11598 


42 


COS 



Note. — A summary of the observations performed under the HST program GO/PAR 11696 for the WISP Survey to date. 
Four additional fields were deferred to a later paper, because they require additional reductions. Note that some targets have 
more than 5 orbits of observations. These are not necessarily all full orbits, however, which is why, for example, a three-orbit 
observation may have more than half the total integration time of a six-orbit observation, and targets with the same number of 
orbits have varying total integration times. 



of "generic" parallel visits has improved the efficiency of 
our observing program by about 50% for a fixed allo- 
cation of orbits. We also requested, and were granted, 
an additional 17 single-orbit parallel observations. Al- 
though these were too short for the purposes of other 
programs, when we observe these new parallel fields in 
only the F140W and G141, we are still able to reach 
comparable depths to the G141 images of many of our 
multi-orbit fields (See Table [I]). 

3. DATA REDUCTION 

Slitless spectroscopy offers a considerable number of 
advantages for our WISP program. The high multiplex- 
ing capabilities offered by the WFC3 grisms combined 
with the observation of a significant number of uncorre- 
lated fields offered by the pure parallel strategy, makes 
it perfectly suited for survey purposes. It allows a wide, 
uncontaminated spectral coverage with stable and accu- 
rate calibration. The ability to combine G102 and G141 
grisms, with a well-behaved overlap gives the unique op- 
portunity to find multiple lines per object at faint mag- 
nitudes. Because of the well-behaved photometry, it is 
possible to combine the spectra with broad-band pho- 
tometry to provide comprehensive and precise spectral 
energy distribution (SED). Moreover, the use of the HST 
offers a much lower sky background level compared to 
the ground. However, to fully exploit the potential of 
this tool, one needs to overcome some limitations. In 
crowded fields, object spectra and orders sometimes over- 
lap, leading to confusion in the spectral extraction and to 
uncertainties in the background estimate. Given the high 
latitude of the WISP fields, crowding did not represent 
a problem in our observations. The spectral resolution 
is determined by the size of the target: for point-like 
sources, this corresponds to the instrument Point Spread 
Function (PSF), while for extended sources it is the spa- 
tial extent of the object in the dispersion direction. 



This basic principle of the slitless reduction is the use 
of the direct images to provide the input catalog of ob- 
jects. The object positions, sizes and shapes are then 
used - in conjunction with a precise knowledge of the 
mapping between the direct and the dispersed frames to 
extract photometric and wavelength calibrated 2D and 
ID spectra from the data using a 3D flat-field cube which 
contains the calibration of each detector pixel for each 
wavelength. This process has been integrated with the 
WFC3 pipehne (CAL WF3 version 2.0) with th e provision 
of the aXe software (jKuemmel et al. I l2009a ^ which has 
been developed to handle all of the HST slitless spec- 
troscopic modes. aXe maps not only the + first order 
spectra but also the zeroth, - first, +/- second and higher 
order spectra which allows the quantitative estimation of 
contamination in the case of overlapping objects. 

We processed all the data (including both grism and 
direct imaging exposures) using the CALWF3 pipeline 
(version 2.0). Direct images were corrected for bias, 
dark, flat-field and gain variation using the best refer- 
ence files obtained from the STScI archive. WFC3 per- 
forms multiple non-destructive reads in the slitless spec- 
tra in MULTIACCUM mode with the sampling sequence 
SPARS 100. This allows cosmic ray rejection in a sin- 
gle exposure sequence. The reduction pipeline fits the 
accumulating signal in the the MultiAccum readouts to 
identify the cosmic rays. In the slitless grism data, pixels 
cannot be associated with a unique wavelength. For this 
reason, CALWF3 applies a unity flat so that no pixel- 
to-pixel correction is applied at this stage. The proper 
flat-fielding is then applied during the extraction of the 
spectra with aXe. 

3.1. Direct Images 

Because our pure-parallel observations do not allow 
dithering, the processed images still contain a significant 
number of bad pixels that were not identified during the 



CALWF3 step. In addition, CALWF3 fails to properly 
flag a fraction of cosmic ray hits, when these were impact- 
ing multiple pixels. We identify and flag these pixels in 
the data quality (DQ) masks using custom IDL routines. 
Finally, the images are combined with the IRAF task 
MULTIDRIZZLE (Koekemoer et al. 2002) using stan- 
dard parameters. This step also corrects for bad pixels 
and cosmic rays, and removes the geometric distortion 
using the latest distortion solution (IDCTAB) files. The 
final drizzled images have accurate astrometry and a uni- 
form photometry across the image. 

The position of the spectra in the dispersed im- 
ages is determined by the location of the objects in 
the corresponding direct exposure. We run SExtr actor 
(|Bertin fc Arnoutslll996[ ) on the drizzled direct images, 
creating independent source catalogs for the FllOW 
and F140W images. For the source detection we used 
the RMS weight maps derived during the MULTIDRIZ- 
ZLE step. Every source in the catalogue is required to 
have at least 6 connected pixels (DETECT_MINAREA 
= 6), with values above 3.5cr (DETECT.THRESH = 
3.5<j). Finally, the catalogs generated for the FllOW and 
F140W images are matched and put on the same num- 
bering scheme. The individual detection in both filters 
allow us to keep objects detected in only one band and 
ensures an accurate astrometry between the grism im- 
age and the reference direct one. Because aXe performs 
the spectral extraction on individual (not distortion- 
corrected) grism frames, the combined catalog is pro- 
jected back to these frames to generate an input object 
list (lOL) of coordinates for each grism frame to be used 
in the aXe extraction. Starting from these lOLs, the 
aXe software computes the reference position and spec- 
tral trace for each object in the dispersed images. 

3.2. Grism Images 

Before proceeding to the actual spectral extraction, the 
sky background must be estimated and subtracted from 
each grism frame. It is possible, within aXe, to estimate 
the background locally by interpolating on the sky re- 
gions above and below the spectral trace. However, we 
decided to perform a global sky subtraction on the entire 
image before proceeding to the extraction of the spectra. 
This method proves more efficient for slitless survey data 
(jPirzk al et al."2004!). where spectra may be be contami- 
nated by the trace of nearby objects, and multiple orders. 

In the slitless mode, each spatial element of the detec- 
tor receives the full sky and telescope foreground emis- 
sion integrated over the grism bandpass. The grisms dis- 
perse the foreground light, but the lack of a slit results 
in a mapping of spatial to spectral domains that rein- 
tegrates uniform sources (e.g., the sky) effectively into 
undispersed emission. A portion of the detector does 
not receive the full zero-order sky and thus there is a 
ramp in the sky in the first one hundred or so columns 
on the detector. The sky foreground in the G141 and 
G102 grisms is a combination of zodiacal light, earth limb 
and thermal emission from the telescope. The processed 
two-dimensional grism images contain a significant signal 
from these sources that must be removed for proper pho- 
tometric calibration, as well as for optimal identification 
of spectral features. 

We constructed a master sky template for each grism 
by making a stack of all of our grism exposures. The in- 



dividual images were first scaled using the image mode, 
then combined with a median stack. The resulting mas- 
ter sky contains low- level residuals from bright sources, 
but these have little impact on our results. Our sky tem- 
plate is being improved as additional independent fields 
are observed and added to the master sky. 

For each grism observation we scaled the appropriate 
master sky and subtracted this from the post-pipeline 
two-dimensional image. We experimented with scaling 
the sky by the mode of the pixel distributions and by 
minimizing the sigma in the sky subtracted image. These 
two approaches gave nearly identical results, and we pro- 
ceeded with the skies scaled using the mode. The scale 
factors varied by more than a factor of two from one field 
to another, reflecting the changing zodiacal light back- 
ground and contribution from the Earth's limb. 

As in the direct images, the CALWF3 processed data 
still showed residual bad pixels and cosmic rays. We con- 
structed a bad-pixel mask in which we flagged all bad 
pixels, including those with low-QE and with high value 
of the dark current. We then interpolated over these 
flagged pixels along the dispersion direction using a sec- 
ond order polynomial. We find that the majority of the 
bad pixels are isolated, and most interpolations were over 
at most 3 pixels. Thus the interpolation step does not 
introduce spurious features into the spectra. Moreover, 
because the position of each bad pixel is recorded in the 
masks, we have checked that the identified emission lines 
(see below) were not impacted by the interpolation pro- 
cess. Finally, we created the final dispersed images by 
averaging together the corrected frames, weighting each 
image by its exposure time. These frames were used in 
our visual inspection. 

We also run MULTIDRIZZLE on the individual grism 
frames. This has the advantage of performing cosmic ray 
rejection and flagging the affected pixels in these frames. 
However the drizzled grism image is not used for extrac- 
tion. The aXe procedure extracts the spectra from the 
individual grism images and all flagged pixels will be ig- 
nored during this process. Another advantage of apply- 
ing MULTIDRIZZLE on both direct and grism frames is 
to avoid potential differences in the WCS data for the 
two images that could lead to a mismatch in the source 
positions in the direct and grism individual images dur- 
ing the extraction. 

3.3. aXe Extraction 

The spectral extraction is performed with the slit- 
less extraction package aXe (version 2.0) developed at 
Space Telescope European Coordinating Facility (ST- 
ECF). The input image list is prepared with each in- 
dividual grism image paired with the direct image which 
is closest in time and its SExtractor catalog. We used 
the latest in-orbit configuration files (jKuntschner et aJTI 
2009a. b") and thermal vacuum (TV) 3 sensitivity cali- 
brations. Some of these reference files will be improved 
for later reductions. We compared several apertures and 
orientation parameter setups to optimize the extraction. 
We used an aperture width of twice the object size listed 
in the input object list, and an extraction direction per- 
pendicular to the spectral trace (see aXe user manual-^^ 
for details). 

^^ http://www.stecf.org/software/slitless_software/axe/ 
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Fig. 1. — Top panel shows a comparison of an individual G141 grism image before and after the reduction process. The raw image 
on the left shows a large number of bad pixels and cosmetic artifacts that are no longer visible on the clean image. Bottom panel: An 
example of different spectral features in our grism images. Red and green rectangles indicate examples of first and second order spectra 
respectively. The blue circles mark the zero order positions, indicating in some cases a contamination to the spectra. Some emission lines 
are also indicated. 

Two-dimensional cutout spectra of each object in 
each exposure are extracted and then combined (with 
aXeDRIZZLE) to create a final high signal-to-noise 
(S/N) ratio spectrum. The combined 2D spectrum is 
then used to extract the final object spectrum. We used 
optimal extraction to increase the S/N ratio of the ex- 
tracted spectra. The spectral contamination by spectra 
of nearby objects is estimated by aXe using the object po- 
sitions and magnitudes in the direct image. We adopted 
the Gaussian emission model, using SExtractor size pa- 
rameters, to approximate the object morphology. The 
simulated spectrum is then extracted and used to esti- 
mate the contamination to each object. 

3.4. Emission line Objects 

Each stacked two-dimensional grism frame was 
searched visually by two or more of the authors. We used 
the dispersion solution from the WFC3 configuration files 
to determine the positions of the aperture beams in both 
grisms. The 1st order of each object is then marked 
with the corresponding aperture beam and the original 
ID number in the lOL. The actual search area is about 
3.3 arcmin^, smaller than the IR channel field of view. 
Because of the shift between direct and grism images, 
spectra on the left part of the field do not have detec- 
tions in the direct image, and zeroth orders cannot be 



identified at the far-right side of the frame. Some exam- 
ples of 1st and 2nd order aperture beams are shown in 
Figure [1] by the long rectangular boxes. 

Besides the overlap between the objects, the most im- 
portant spurious feature is the zero-order image. Because 
the light is not dispersed, the zero-order image could be 
easily confused with an emission line. We derived the po- 
sition of all the zero-orders in the grism exposures, using 
the positions of all the sources in the direct images. We 
marked the zero-orders in the grism exposures to be able 
to exclude them during the visual inspection of the grism 
images. In Figure [1] we show an example of this proce- 
dure: the blue circles show the zero order location of 
some of the objects detected in the direct images. With 
all possible contaminants identified in the grism images, 
we proceeded to compile a list of all emission-line source 
in each field. We compared the G102 and G141 spec- 
tra of each field side-by-side, to aid in the identification 
of genuine emission-lines. In many cases the identifica- 
tion of the lines is unambiguous. For exampe, galaxies 
with Ha often also have detectable [SII] 671 7/6731 and 
sometimes [SIII] 9069,9532 emission (Figure EJ. Strong 
[OIII] 5007,4959 emission in the G102 grism is easily rec- 
ognized on the basis of Ha in the G141 and often H/3 
in the G102 spectra, even though the two [OIII] lines 
are blended (see Figure [3]). At higher redshifts, strong 
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[OIII] 5007,4959 in the G141 grism is confirmed by faint 
H/3 emission. In some cases [Oil] 3727 can be seen in the 
G102 spectra for sources with strong [OIII] 5007 in the 
G141 spectra (Figure |4]). 

There are a large number of objects with only a single 
faint emission-line in the full spectral range covered by 
the two grisms. It it most probable that these are weak 
Ha emission lines for which we are unable to detect other 
lines. However, since we see many [OIII] /Ha pairs up to 
z ^ 1.5 we expect, in absence of other indication, to have 
a fraction of the single lines in the G141 that would cor- 
respond to [OIII] whereas Ha has moved out of the G141 
grism. The [OH] 3727 line is seldom recognized as a red- 
shift indicator by itself, since in this case it is more likely 
to be an Ha line at low-z rather than a high-z [O II] line. 
There are some objects, however, for which we believe 
the single line in the G141 spectrum is [OIII] 5007,4959. 
This is primarily based on the lack of continuum in the 
G102 spectra, indicative of the fall off in the continuum 
below the Balmer and 4000 angstrom breaks. 

Each emission feature was assigned a confidence flag, 
ranging from to 4 in order of decreasing confidence. 
The highest confidence level (class 0) was assigned to 
objects with multiple emission features and unambigu- 
ous redshifts. Confidence classes 1, 2 and 3 were as- 
signed to single features with unambiguous identifica- 
tions (class 1), good quality (class 2), moderate quality 
(class 3), and uncertain (class 4, which are not included 
in the present analysis). Similar confidence classes have 
been use d in gro und-based redshift survey s (e.g;. CFRS, 
lie Fevre et al. I ([2005); GDDS, Abraham e t al. I (|2004)) 
and have been found to be useful in restricting various 
types of analysis to high confidence objects. Monte Carlo 
simulations will be used to assess the flux completeness 
and reassess the confidence levels when the data set is 
complete. 

The final object catalog contains a measurement of 
each emission line, performed using a custom IDL fitting 
tool, and relevant information from the object extraction. 
The line fitting uses the pixel position of the identified 
lines, applies the grism dispersion solution to find the cor- 
responding wavelength and therefore the redshift of the 
object. The fitting model consists of a polynomial con- 
tinuum and a gaussian (a multi-gaussian model is used to 
deblend the Ha-[Sll] and [Olll]A5007,4959 A lines). The 
program fits the identified lines by least squares, and tries 
to fit potential lines that could be present in the spec- 
trum at the correct redshift with a lower significance of 
typically 2 <j. The entire process is visually monitored 
by displaying the line, the fit model and the robustness 
of the fit. The fitting flux uncertainties are derived from 
our aXe extraction and calibration. As an independent 
check, a large subset of these same emission lines were 
also measured by hand, using the SPLOT task in IRAF. 
The resulting splot measurements generally agree with 
the automated measures to within the uncertainties. 

We investigated the accuracy of the grism data cali- 
bration. First, we compared the fluxes of lines detected 
simultaneously in both G102 and G141 grisms, in a wave- 
length domain restricted to a throughput higher than 
10%. In most cases the flux agrees to better than 10% 
between the two grisms. Second, we compared the agree- 
ment between wavelength calibrations by using the over- 



lapping region around 1.15 micron, and the redshift de- 
termined independently from different lines detected in 
G102 and G141. We observed a typical offset of 2 x 10"^ 
microns that is within the uncertainties of the redshift 
determination in the individual grisms. 

3.5. Grism Sensitivities 

Due to the pure parallel nature of our program, we 
have no control over the location and timing of our ob- 
servations. Because of this, we are subject to a variety of 
visit durations (typically 3-5 orbits), number of readouts, 
and backgrounds (zodiacal light and earthshine), all of 
which affect the sensitivity of our observations. In or- 
der to determine the depth of our spectroscopy, we have 
measured the background rms (in electrons per second) 
of every final reduced grism image. This rms is then di- 
vided by the sensitivity (as a function of wavelength) and 
divided again by the dispersion (in A pixel" ^) to give the 
la noise as a function of wavelength. In Figure [5] we plot 
the 5<j continuum sensitivity in a five-pixel wide aper- 
ture (in the spatial direction) for each pointing. We used 
three pixels width in the dispersion axis to estimate the 
line sensitivities. Each field has been scaled by (exptime) 
^/^ so that the relative depths are comparable. Though 
most observations agree with the predicted depths from 
the WFC3 exposure time calculator, there is a large dis- 
persion in the depths due to background variations, and 
any given field can vary substantially from the mean by 
typically 20% (up to ~ 50%). We also plot in Figure 
[6] our line sensitivity as a function of depth. In both 
grisms, we reach a median depth of about ~ 5 x 10~^^ 
ergs s~^ cm~^. 

In summary, for the typical exposure times of the sur- 
vey (- 6400 s in G102, - 2700 s in G141) we are able to 
detect (at 5cr) compact lines with / = 5x 10~^^ ergs s~^ 
cm-2 at 0.95 < A < 1.17 A and 1.3 < A < 1.65A. 

4. EMISSION LINE RESULTS 

To date, we have selected emission lines over 19 WFC3 
fields, and measured redshifts and line fluxes from the 
highest confidence emitters. While incomplete at the 
faintest lines, we clearly reach emission line fluxes down 
to a few xlO~^^ erg s~^ cm~^, often in the presence of 
a second brighter line. At the highest confidence levels, 
our preliminary analysis yields 328 (S/N > 5) emission 
lines from 229 objects in our 19 WFC3 fields. Simula- 
tions will be required to determine our incompleteness 
due to extended emission, confusion, both high and low 
equivalent widths, while accounting for the non-uniform 
survey sensitivity shown in Figure [5l We list the iden- 
tifications of these lines in Table 2. We note that the 
fraction of sources with multiple emission lines increases 
if we include lines detected at 3<j to 5cr. 

We find 185 galaxies with likely or definite Ha emis- 
sion. Amongst these, 56 have secure identifications 
based on multiple emission lines. We also identify 
single-line objects as Ha emitters, as was the case 
with NICMOS g r ism spectroscopy (iMcCarthv et al. 1 
119991 : iHicks et al.l 120021 : IShim et al.l l2009f ). This as- 
sumption has generally worked in ground-based narrow- 
band IR detections of singjle emission line s (e.g; . 



Bunker et al.l ll995MMalkan Teulitz fcMcLeanI 119961 : 



van der Werf et al.ll2000l : ISobral et al.ll2009f ). 
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Fig. 2. — Examples I: low-redshift emission line objects at 0.4 < z < 0.8. For each object we show the direct image cut-out (5^^ x 5^^), 
the 2D G102 and G141 grism spectra and the ID extracted spectrum at observed wavelength. In this category, Ha falls in the G102 grism. 
We can see several emission hues: [S II]AA6716+6732 A, [S III]A9069 A, [S III]A9532 A and He I A10830 A. 
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Fig. 3. — Examples II: emission line objects in the redshift range 0.6 < z < 1.3. For each object we show the direct image cut-out 
(5^^ X 5^^), the 2D G102 and G141 grism spectra and the ID extracted spectrum at observed wavelength. This category includes mainly 
Ha[0 III] hue emitters, but we also observe [O II], [Ne III], H/3, H7 and [S II]AA6716+6732 A. 
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Fig. 4. — Examples III: high-redshift emission line objects between 1.3 < z < 2.3. For each object we show the direct image cut-out 
(5^^ X 5^^), the 2D G102 and G141 grism spectra and the ID extracted spectrum at observed wavelength. At these redshifts, both Ho; and 
[O III] are shifted to the G141 grism. In many cases we also observe the [O II] emission line in the G102 grism. 
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TABLE 2 

Emission lines in the WISP Survey 



1.0 1.2 1.4 

Wovelength [/i.m] 

Fig. 5. — The 5a continuum (top) and line (bottom) sensitivities 
as a function of wavelength. The depth of each field is plotted in 
gray and scaled (by exptime-*^/^) to the "nominal" survey exposure 
times. The thick solid lines denote the median field depth. The 
blue diamonds are the values given by the exposure time calculator 
assuming "average" earthshine and zodiacal light backgrounds. A 
five-pixel (0.65^^) wide aperture was used for these estimates and 
the line sensitivities assumed three pixel widths in the dispersion 
direction. More compact sources will give slightly better sensitivi- 
ties. There is significant scatter from the median due to variations 
in background levels. 
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Fig. 6. — Cumulative area as a function of depth. Each field has 
an effective area of 3.3 arcmin^ which is added to the cumulative 
area at the value of its five sigma line sensitivity (in a 5x3 pixel 
region). The depth varies as a function of wavelength so we choose 
to use the depths at the central wavelengths (1.0 /im for G102 and 
1.4 /im for G141). The depths are better at longer wavelengths and 
worse at lower wavelengths. The median depths in both grisms are 
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Note. — The main emission lines 
with a signal to noise ratio S/N > 5 
detected in the WISP Survey. 

Other possible line identifications would almost al- 
ways result in a brighter line that we would also de- 
tect. The most significant possibility of a single emis- 
sion line that is not Ha, would be [OIII] 5007+4959. 
This doe s occasionally happen in ground-based near-IR 
searches (|Teplitz et al., |[l999f ). Not surprisingly, for the 
remaining 56 objects with multiple bright emission lines, 
the most common combination is Ha + [OIII] 4959/5007. 
The redshift distributions for the total sample is shown 
in the left panel of Figure [71 while the right-hand panel 
shows the redshift distribution for the Ha, [OIII], and 
[OH] emitters separately. 

We next compare these first results to our previous 
NICM OS Ha survey (McCarthv d~^[T999HShim et al.l 
l2009f ). With WFC3, we see Ha emission from 0.3 < z < 
1.5, whereas with the NICMOS G141 survey we were sen- 
sitive to Ha from 0.7<z<1.9. In the left panel of Fig- 
ure [TOl we plot raw number counts for the Ha and single 
line emitters separately (red and blue points respectively) 
then for the total (Ha + single-line) sample (black). We 
then compare total counts from both surveys, in the over- 
lapping redshift range of 0.7 < z < 1.5 (right panel). At 
the brighter end, we find a good agreement between the 
WISP survey and the NICMOS grism results, even with- 
out including the many lines which we have identified 
with a lower signal-to-noise ratio. But WFC3 starts to 
detect more galaxies than NICMOS at / ~ 3 x 10~^^ ergs 
s~^ cm~^. Towards the faint-end and the completeness 
limit of NICMOS, the surface density of Ha emitters de- 
tected with WFC3 is a factor of 3 higher. This is because 
our survey is more complete for extended line-emission 
and our spectra are sensitive to a very wide range of emis- 
sion line equivalent widths, from cases where the contin- 
uum is hardly detected (EWo6s > 300 A, see Figure [9]), 
to low contrast features with EW as low as 15 A for lines 
that satisfy S/N = 5<j. 

As Figure [8] illustrates, we routinely detect emission- 
lines in galaxies with luminosities well below L*. Many 
of these very faint galaxies have strong [OIII] 5007 emis- 
sion and relatively weak [OH] and H/3 emission, sugges- 
tive of low metal abundances. The near-IR photome- 
try provided by our FllOW and F140W images, and 
ground-based follow-up (where practical), will allow us 
to extend pioneering w ork on the mass-m et allicity rela- 
tion at 0.5 < ^ < I (Savaglio et al."2005) and 1 < z <3 
(|Erb et al. I l2006f ) to significantly lower stellar masses. 
The mass-metalicity relation is a critical diagnostic of 
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Fig. 7. — Le/t - The redshift distribution of our total sample of emission-line galaxies in 63 arcmin^ of the WISP survey. Filled histogram 
represents our secure redshift determination (class and 1), whereas on the top, the hatched part shows the uncertain single-emission-line 
objects assumed to be Ho; emitters. Right- The redshift distribution of our galaxies according to their emission lines. Again the hatched 
part of the Ha emitters represent uncertain redshifts. 

ies is -2 0.75, one and ha l f mag;nitudes below L"^ for z 
-1.0 (Chen et all [20031 : iColbert et al.l[2006h . At the 
lower redshifts, especially for compact objects, our sur- 
vey is able to detect lines in extraordinarily faint galaxies 
{Mr > —19), and star formation rates less than 1 solar- 
masse per year. 

4.1. Absorption Line Spectra 

Although we have so far limited the analysis to emis- 
sion lines, the broad wavelength coverage and high sensi- 
tivity of the WISP survey leads also to a wealth of infor- 
mation from continuum spectra. The 4000A break and 
the convergence of the Balmer series at 3650A, provide 
two of the best low-resolution diagnostics for redshift and 
stellar populations. Continuum spectral diagnostics, par- 
ticularly at z > 1, provide constraints on the mean age of 
the stellar population and the time since the last major 
episode of star formation. Spectroscopy of red galax- 
ies from the ground in the visi ble is limited to objects 
with z < 1.8 and / < 25 (e.g. 'McCart hv et al. I l2004l : 
ICimatti et alll2004l : [Pohertv et al. 2005). Near-IR spec- 
troscopy can prob e to hig;her redshifts, but this is quite 
challenging (e.g. Kriek et ahl 12008.) and has been lin i- 
ited to small samples to date (e.g. iDohertv et aT]|2006f ). 
The WFC3 grisms are providing low-resolution contin- 
uum spectroscopy of red galaxies in large numbers with 
great sensitivity (e.g. Ivan Dokkum fc Brammer I l2010l ). 
In Figure [12] we show an example of a post-starburst 
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Fig. 8. — H_140-band magnitude as a function of redshift. Secure 
redshifts are shown as black circles and single line (assumed to be 
Ho;) emitters as red crosses. We also show the apparent H_140-band 
magnitude for L"^ (and 0.1 L"*") galaxies derived from rest-frame 
near-IR (at low- 2;) and optical (at z > 0.7) luminosity functions 
(Bell et al. 2003, Faber et al. 2007, Shapley et al. 2001). 

galaxy formation and evolution models as it is sensitive 
to a number of physical process (e.g. infall, SNe feed- 
back, IMF). 

In Figure [TT] we plot on the left panel the observed 
star formation rate distribution of o ur Ra sample de - 
rived from the Ha lu minosity using t he iKennicuttl ll998l) 
conversion. Following IShaplev et al.l (|2Q05). we corrected 
the Ra fluxes for 20% average [Nil] 6584+6548 contribu- 
tion, but no extinction correction is applied. On the right 
panel we plot SFR(HQf) versus rest-frame R-band lumi- 
nosity. For rest wavelength R magnitudes we used the 
FllOW photometry if the observed wavelength of Ha is 
< 1.25 /im {z < 0.9) and the F140W photometry for Ha 
> 1.25 /im {z > 0.9), applying only the {l-\-z) correction 
(without any k-correction). There is no extinction cor- 
rection for the R-band magnitudes. The median SFR is 
4 M0/year, and the median M^^ of line-detected galax- 



galaxy at redshift z = 1.68. A clear Balmer break is seen 
in the G102 spectrum, as well as absorption lines in the 
G141 part. A stellar population model is overplotted, 
giving an age of 600 Myr and an extinction of Ay = 1.2. 
This object has been pre-identifled searching for galax- 
ies with H_UOab < 22 and J_110 - H_UO > 0.6. This 
selection isolates bright passive galaxies at z > 1. 

5. DISCUSSION AND SUMMARY 

We have presented the flrst analysis of the WISP Sur- 
vey of slitless spectroscopy with the NIR channel of 
WFC3, taken in pure parallel mode. In the flrst 19 flelds, 
we surveyed 63 arcmin^. At the highest confldence level 
we have detected 229 emission line objects. We reach 
typical 5cr sensitivities of / ~ 5 x 10~^^ ergs s~^ cm~^ 
for compact lines. 
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Fig. 9. — Examples IV: very high equivalent width emission line objects. A strong emission line is observed in this group of galaxies with 
an extremely faint continuum which barely visible in the 2D grism images. For each object we show the direct image cut-out (5^^ x 5^^), 
the 2D G102 and G141 grism spectra and the ID extracted spectrum at observed wavelength. 



13 



§" 1.000 






0.100 



0.010 



0.001 



t 


' ' 1 

t 


i 




1 . 

Ha (Confirmed) • : 

Single Line • ■ 

Ha + Single Line • . 




i 




I 


■ 


r 




{ 


\ 


, 


- 










' 


' 












' 










! 


' ~. 




. . 1 








■ 







' ' 1 




1 


1 ■ ■ ■ 1 . 

WFC3. : 

MICMOSn ■ 


£;:^ 


* 








■ 


8" 1.000 


- 


* 






- 


arcmin /d( 

p 

o 
o 


r 


[] 


f 


* 


- 


E 0.010 








[] 

( I 


=5 










f- 


0.001 




. . 1 








. 



Line flux (1 0'^^ ergs s'^ cm'^ 



10 



1 10 

Line flux (10'^^ ergs s'^ cm"^ ) 



Fig. 10. — Left -Number counts of Ho; emitters in the WISP survey. Confirmed Ho; emitters are shown in red points, single-fine emitters 
in blue and the total counts in black. Right- Number counts of Ha and single line emitters in the WISP survey (black points) compared 
to the NICMOS Parallel grism survey (green squares), tabulated in the overlapping redshift range of 0.7 < z < 1.5. Counts are given per 
d{logf) = 0.25. 
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Fig. 11. — Lejt - Observed star formatio n rate dis t ribut ion of our Ho; emitters. The Ho; fiux is reduced to account for 20% of [N II] 
contamination and converted to SFR using Kennicutt (1998) relation. The sample is split into high-z (shaded area) and low-z (hatched) 
sub-samples, where the separation between the two categories is z=0.675. Right - SFR as a function of absolute R-band magnitude derived 
from FllOW and F140W photometry for Ho; emission below and above 1.25 microns respectively. Black points represent confirmed Ha 
emitters and red crosses the single-line galaxies. 
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Fig. 12. — Example V: Spectrum of a post-starburst galaxy at z = 1.68. A clear Balmer break can be seen in the G102 spectrum between 
1.0 < A < 1.1/im as well as Balmer absorption lines in the G141 spectrum at 1.17 and 1.3 /im. Plotted (in red) is a Bruzual &; Chariot 
(2003) stellar population model with an age of 600 Myr (e- folding time of 100 Myr), Ay = 1.2 mags (assuming a Calzetti reddening curve), 
and a Salpeter IMF. 

their first major star formation episode. 

The sensitivity and resolution of the G102 grism also 
provide the opportunity to search for Lya emission lines 
in addition to the rest-frame optical emission lines dis- 
cussed here. We have excluded such analysis from the 
present paper given the limitations of the initial, vi- 
sual selection of faint emission lines. After further anal- 
ysis, WFC3 will offer a significant new approach to 
high-redshift {z > 6) searches for Lya. Though exten- 
sive ground-based narrowband searches for Lya-emitters 
have been performed up to z = 7, with the most dis- 
tant confirmed source at 2; = 6.96 (|Ive et al. I I2QQ6). it 
is difficult beyond this redshift to find wavelengths with 
both high sky transmission and low sky emission. Thus, 
many surveys a re limited to a very small range at z = 7.7 
OT z = 8.8 (e.g. ISobral et al. 2009), without any reliable 
Lyof emission-line candidates. 

In addition to the WISP survey, two large observing 
programs with the WFC3 grisms will contribute sig- 
nificant results. First, a G141 survey of the GOODS 
North fields is underway (GO-11600; PI=Weiner), reach- 
ing depths comparable or slightly better than WISP. 
Secondly, the CANDELS multi-cycle treasury program 
(GO-12060 through 12064; PI=Faber and Ferguson) wih 
obtain deep spectroscopy of selected pointings targetted 
to follow-up supernovae found in direct imaging. The 
WISP survey has some limitations compared to these 
programs due to the parallel mode, including the lack of 
dithering, but our first analysis demonstrates the great 
efficiency of this use of the telescope. WFC3 parallels 
can obtain large spectroscopic samples of emission-line 
objects during the peak epoch of star formation, free 
from the most common selection biases. 

We gratefully acknowledge the dedicated efforts of sev- 
eral staff members at the Space Telescope Institute to 
make the new Grism Parallel observations as successful 



While wide-field narrow-band imag i ng surveys from 
the ground (e.g. ISobral et al.l 120091 : iMomcheva et al.l 
I2OIOI ) yield samples of Ha emitters that are larger than 
ours and near-IR slit spectroscopy of selected targets can 
yield emission-line ratios to faint levels (e.g. Pettini et 
al. 2001), the WFC3 grisms offer a unique opportunity 
to obtain continuous spectral coverage in the 0.8 — 1.7/im 
region with high sensitivity. Our early results show that 
we can identify star forming galaxies with space densities 
comparable to, or larger than, those identified with color 
selection techniques (e.g. BM/BX or LBG selection; 
Adelberger et al. 2004) to similar star formation rate 
thresholds (Shapley et al. 2005, Erb et al. 2006). The 
typical emission-line galaxy at z > 1.3 in our sample has 
strong [OIII] 5007,4959 emission and a large [O IIl]/H/3 
ratio, similar to the ratios seen in the small number of 
LBGs at z > 1.5 with rest-frame optical spectra (Pettini 
et al. 2001, Teplitz et al. 2000, Lemoine-Busserolle et al. 
2003, Maiolino et al. 2008, Mainline et al. 2009). Cou- 
pled with the high [O III] /Ha ratio and weak [O II] 3727 
these suggest low metal abundances. Future papers from 
our survey will examine line ratios and their implications 
for abundances and ionization sources. With the support 
of our ongoing ground-based, follow-up observing pro- 
grams, we will be able to measure the mass-metallicity 
relation at crucial intermediate redshift s, to derive the 
Ha luminosity function at different redshifts and to study 
the evolution of the star formation rate density. 

Many of our fields contain emission-line objects with 
extremely high equivalent widths (cf. Figure [9]). These 
strongly star forming emission-line dominated objects 
are reminiscent of the ultra-strong emissio n-line ( USEL) 
galaxies at < z < 1 disc overed by Kakazu et al. I (J2007) 
(see also iHu et al71 l2009f ) and the "green pea" galaxies 
found in SDSS by Cardamone et al. (2009). These ob- 
jects have metallicities comparable to the lowest observed 
in local sources and may be dwarf galaxies undergoing 
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